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SUMMARY 


% 


The theoretical predictions of several versions of the NASA- Lockheed multielement 
airfoil computer program are evaluated. The computed results are compared with 
experimental high lift data of general aviation airfoils with a single trailing edge flap, 
and of airfoils with a leading edge flap and double slotted trailing edge flaps. Detailed 
test-theory comparisons are reported for the GA(W)-1 with a 30% chord flap and a 
Boeing four-element high-lift airfoil. Theoretical and experimental data comprise global 
airfoil parameters such as lift, pitching moment, and profile drag, as well as surface- 
pressure distributions, boundary layer integral parameters, skin friction coefficients, 
and velocity profiles. The evaluated versions of the computer program include the 
baseline version that was available from NASA in June 1976, and three versions that 
were developed by The Boeing Company in a joint NASA-Boeing effort. The results of 
this evaluation show that the contract objectives of improving program ^reliability and) 
accuracy have beemmetwft?.*- 
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This document reports on an evaluation of the predictions of several versions of the ) 
NASA-Lockheed computer program for two-dimensional multielement airfoils. The f 
original version was developed by Goradia and his coworkers at Lockheed Georgia 
under the sponsorship of the NASA-Langley Research Center (ref. 1). The program was 
later modified extensively in order to improve its predictions for the different types of 
high-lift airfoils. Many improvements, mainly in the area of the potential flow 
calculation, were mad e by researchers atjthe Langley 'Research Center (ref. 2) . Recently ; { d 
members of the Aerodynamic Research Group of the Boeing Commercial Airplghe| 
Company supported by numerical analysts of Boeing Computer Services ;(Bl| s)i 
substantially modified the aerodynamic model, its numerical implementation, and^ti 
computer code. , 


This work, reported 4 in two engineering documents (refs. 3 and 4) and in this 
supplemental document, was sponsored jointly by NASA and ! The Boeing Company. In 
particular, the evaluation of the computer program by comparison with recent 
experimental high-lift data was a joint NASA-Boeing effort. 

The objective of the evaluation was to learn as much as possible about the range of 
applicability and the deficiencies of the aerodynamic model. For that reason those sets 


of experimental data were preferred and investigated in detail which not only offered 
global aerodynamic force and moment coefficients but also information on pressure^, 
distributions and boundary layer parameters. i 


It should be emphasized that the test theory comparison reported in this document is by 
no means complete or exhaustive. However, enough information has been accumulated 
during the evaluation period to justify the publication of a separate document. Revisions 

of the document will Be, made later as .new .versions of the computer program and 

* .»■ . * ~t. -«• » ,-s. . * .» *• .4? '.ft "V *. s§( *■ * 


additional experimental dkta become available.* 
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PROGRAM VERSIONS | 

The following defines the various versions of the computer program emphasizing their . 
differences. J 

... • ■ 

BASELINE VERSION ' | 

This version was supplied to Boeing by NASA-Langley prior to the beginning of the \ 
contract work. A discription is not available but most aspects of its underlying 
aerodynamic theory are discussed in references 1 and 2. The main assumptions of the 
aerodynamic model are: [ 


• Flow is two-dimensional and subcritical 

• Flow is attached to the airfoil’s surface 

* t *. 

• High-lift airfoil consists of up to four components 

• Geometry of the high-lift airfoil is discretized by up to 165 surface points 
VERSION A 


' " h \ 


This is the baseline version of the code made operational for airfoils with negative i 

' £ 

overlap of neighboring components. 


VERSION B 


Sf o 


'T».s version is thoroughly documented in reference 3. The following areas differentiate 
Version B from the baseline version. • , 




The prediction method of Nash for ordinary turbulent boundary layer flow is 
replaced by the method of Nash and Hicks |< ref. 5). 

The drag prediction method of Squire and Young href, fit replaces “the previous-^. 

' v (■ jfl * ' * 

pressure and skin friction, integration. > * * 




* 


Several logical errors in the baseline versionja^corrected. 


Jf 


% § * 


VERSION C 


^ 8?s- 

. 4 # * 


#' * * 

0 <? * * 


1 1 


r # ■* ■ 


■*» 


The following discribes the differences between the aerodynamic theory and the Version h 
B theory of the NASA-Lockheed program: (ref.. 4).. ' * * * * ” . ,j v‘ * * Ti 


• | The method used to represent the displacement effect of viscous layers is replaced jj|| 

I by the surface transpiration method. This uses an equivalent •distributiojb ;,bf'\ * 

| sources along the airfoil surface and along the wake centerlinesyh 


^ y- - 





" "W a *-”' y ~ -♦» aty ■■ , l uy — »» « m- 1 * *<»** -'**' «" -ii -™«w 

The flow model of the region containing a potential core is modified. Wake 
properties are calculated using the lag entrainment method of Green (ref. 7) and 
transition from laminar to turbulent' boundary layer flow may take place in the 

, . /- . ■ f>. ... * i* 0 fe -> 

*core region. „ 

<R •?) "# " ♦ « t* * ’• 

By using a modified version of Goradia’a confluent boundary layer model an 
attempt is made to predict’ the onset *of ft separation ‘of confluent boundary layers. 
The modification utilizes 5 Coles’ fwo-parameter velocity profile for the wall layer 
(ref. 8). *• * . * •• ' « • * - '■ * * ** *• " * * v 


The high-lift airfoil may consist of up to 10 airfoil components; 3 ^ * ***' ' 

N t (* , *#• 4-t 0' ^ *r C \ " 

The computer code of this version uses a dynamic storage allocation thereby 
removing the limitation on the number of geometry points representing 'airfoil 


surfaces. 


5£-' 


.. if 

TEST CASES, 


*Ss- 


' ' ;'#! 

. 3-? 

■# / 


*-* 

' v$ ‘ ' if- -#>"■ .St 
if •* .* ■ *: 


The geometries of all airfoil configurations that were used for the evaluation are shown 
in figure 1. Geometric details of these airfoils such as flap settings, gap, and overlap of 
airfoil components are listed in tables 1, 2, and 3. Furthermore these tables contain the 


l ight condition, i.e., angle of attack, Mach number, and Reynolds number, at which the 


Lest theory comparisons were conducted. Figure 2 illustrates the definitions of thejnost 


important geometry parameters. 

• ' — .wSU:.' * • ' . . .im M fWjl 
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Table 1. — GA (W)-1 Airfoil Cases 


(a) GA (W)-1 Basic Airfoil, c re f -1.917 ft 




| a = -4.11°,j0°, 4.17°, 8.02°, 12.04° 

Moo = 0. 1 5 “ 

R n = 6.3 x 10 6 


(b) GA(W)-1 With 30% Chord Trailing-Edge Flap, c ref = 2 ift 



■H 

j Gap 
j % chord 

; Overlap 
; % chord 

'& 

\ degree 

; io 

: 2.5 

j 7- 1 


-4.0, 0.2, 5.2, 10.3, 12.8 

20 

2.5 

1 4.7 


0.2, 5.2, 10.3, 12.8 

30 

P 2.5 

2.1 


-4.0, 0.1, 5.2, 10.3, 12.8 

40 

2.6 

-0.85 


2.7 


! M oo = 0.13 
: R n = 2.2 x 10 6 











\ 
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I Table 2. — Additional Two -Element A irfoil Cases 


(a) Foster's Two-Element Airfoil, c re f = 3 ft 



i Gap 2.5% chord 
; Overlap 4.3% chord 


i 5 F = 30° 

! M m 1 = 0.18 (?) 

; R n = 3.6 x 10® (?) 

i d = 0° 


\ 

(b) Ljungstrom's Two-Element Airfoil, c re f = 2.133 ft 



f Configuration 1 , 32% chord flap 


§F 

= 20° 

Gap 2% chord 

Moo 

j= 0.12 

! Overlap 1.54% chord 

Rn 

= 1.83 xIO 6 


a 

= 0° 
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Table 3. — Four-Element Airfoils 




Slat 

j Flap 1 j 

| Flap 2 

!. 5 S j 

| Gap 
f %c ref 

...... 

Chord 

% c ref 

: 1 S F 

Gap 

| %c ref i 

j Chord 
j %c ref 

fi s F ~1 

; Gap 
'■I % c ref f 

Chord 1 
%e re f I 

: 53 ° 

! 1.4 .j 

• " . ;i 

11' 11.2 | 

1 11° : 

: 4 i 

i 2.8 

j 19.4 

| i 13° ! 

| 0.5 I 

j 9.9 


\ Moo = 0.16 

| R|M = 2 * 106 

a = 4.4°, 8.4°, -12.45°, 16.6°, 20.6° 


J (b) Ljungstrom Four-Element Airfoil, C re f = 2.133 ft 
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boundary layer or wake momentum thickness 
boundary layer, thickness on upper surface 




f w * >. ; . ^abbreviations ' - 


FT , 


»LS. 




? T, 


' " -%?' / ' 


fixed transition 

% 

laminar short bubble 
1 separation ■„* ,- 
... : transition. ^ 


.K 





10 







EVALUATION OF PROGRAM PREDICTIONS 


vi 


The comparisons of theoretical and experimental high-lift data are discussed in the j 
same order in which the analyzed airfoil configurations are .listed in figure 1 and s 
tables 1, 2, and 3.' 1 *' * 


M 


SINGLE AIRFOIL 


§ jt € £ 1? 

xki 'M p & '*+ 


BASIC GA(W)-1 AIRFOIL 


The basic GA(W)-1 airfoil was chosen - to test the program capability of predicting 
performance characteristics of single airfoils. Figures ! 3 and 4 contain theoretical lift, 
pitching moment, and drag curves and their comparison with experimental data of 
McGhee and Beasley (ref. 9h Both Version A and the new Version C programs predict^J 
identical lift, and moment curves that in turn agree well with measured GA(W)-1 data 
up to the onset of trailing edge stall at about 8° angle of attack. 

Considerable differences between all drag polars are observed in figure 4. Version A, 
utilizing an integration of surface pressure and skin friction in the prediction qfiprofile* 
drag, gives the highest drag coefficients. Version C, applying the Squire and Young,. . 
formula, offers drag values that are lower than the corresponding experimental drag'b 
coefficients. The lack of agreement of die three drag polars emphasizes the fact that • 
even for single airfoils at low speed the problem of obtaining acourate^dra ^ % 
computations is not yet solved. *' ‘ J 




re 


Surface pressures of the GA(W)-1, computed by .program Version C and plotted ini ft; 

5 at two different angles of attack, agree reasonably well with their experimental f 
counterpart. At the lower angle of . attack the flow is attached to the entire airfoil , I, 
surface; whereas, at the higher angle of attack a laminar short bubble with subsequent, 
turbulent reattachment of the boundary layer is indicated near the upper surface.? ! 
leading edge, and turbulent boundary layer separation is predicted theoretically j near/ 
the upper surface trailing edge. The latter prediction is confirmed by the experimental? 
pressure distribution which shows a constant pressure downstream of the theoreticdf$|| 
point of separation. 


•4P ! 

A 


In all figures of this document, the arrows and the symbols S and LS refer toHheoretical 
points of turbulent separation and laminar short bubbles, respectively. The symbol FT 
indicates the experimental trip strip location which is specified as a fixed transition 


point in the computer simulation. The symbol T denotes a computed free transition 
point. \ ■ ■ ‘ ■ 


FLAT PLATE 


•I i * 



1. 


Figure 6 shows a test theory comparison of a different kind that was used to judge the < 
quality of the wake flow calculations of Green’s lag entrainment method (ref. 7). The 3 " 
experimental data of the wake behind a flat plate at zero angle of attack were taken 
from test case 14 of the NASA-Langlev Conference o n,Free T urbul ent Shear^lows (refst. 


li 









O i Expe riment ^ McGhee & Be asley 



, Figure 5. — Surface Pressures of GA (W)- 1 Single Airfoil 




o 

Experiment: McGhee & Beasley 

a = 8.02° 
Moo =0.15 


Version C 

;R n = 6,3 x 10 6 

FT 

Fixed transition 


S 

Separation 


LS 

Laminar short bubble 



Figure 5. — f Concluded f 







10 and 11); ana appear to be the only available test results of two-dimensional wake 
characteristics at low-subsonic speeds. The agreement of theoretical and experimental 
wake parameters seems to be adequate for the present purpose, but better results for 
the wake displacement effect immediately downstream of the trailing edge might be 
r desirable. . 

riinurtiitim' ii 'fotti’n — o&af?. — , 
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TWO-ELEMENT AIRFOILS 
GA(W)-1 WITH 30% CHORD FLAP 


The GA(W)-1 airfoil, with a single 30% chord trailing edge flap, served as the principal 
test case for this type of general aviation high lift airfoil'. The experimental data were 
measured by Wentz, Seetharam, and Fiscko (refs. . 12, 13, and 14). The data' include, 
global airfoil parameters, detailed surface^, pressures,' and boundary layer- 
characteristics. «■- ’ * ?..? •« * *• »-■ >' '# <* * '■ * • .... 

Figure 7 shows -the' lift and pitching moment characteristics of this%irfoil with a flap 
deflection of 10°. The computed data of Version C agye§. veryjwell with the'experimental 
results in the prestall angle of attack range; Version B slightly mispredicts lift^ and 
moment curves. Figure 8 presents the lift curves at three different flap settings up to 
- 30° flap angle,. j ■- -- - " * ? % _ * % v ; i * ■ 




w- 


Figure 9 shows the corresponding pitching moment characteristics,* The ^agreement of 
the theoretical data with the experiment is only satisfactory at lower jangles of attack 
and lower flap angles. The differences' at higher angles of attadk arer due to 
trailing edge stall, a phenomenon thatSis InoV modeled by the computer programs * 

* > • r * 1 l r , » « 4 * a --1 #' + ■* ? ; \ * 

fThe discrepancy between theoretical and ^experimental lift coefficients at the 
30° flap angle and small angles of attack is unexpected in view of the good agreement 
between the theoretical and experimental pressure distribufciohs (fig. 10)., Furthermore,- 
the experimental pressures do not indicate separation on the upper surface of the flap at 
small angles of attack. Hence, trailing, edge stall can not be responsible for the observed 

disagreement? * * * * * * . „ « f & *• . . ,v 

«*-»»****'• + $ « m ■■ • <* * * ■■■ . . 

i Figure il shows a large discrepancy .between theoretical and experimental drag, 
(polars. Presently, one can only speculate about the cause of 'the observed differences. 
None of the drag polars shown in the figure is believed to be correct for the foll owing 
reasons: the theoretical values of lift and drag were computed on the basis of a model 
with many simplifying assumptions; the validity , of gsopie^oi^ these assumptions is 
questionable. A case in point is the computation of profile drag applying thelSquire and 
Young formula, which is a crude approximation for high lift airfoils and, therefore^ is 
not expected to produce very accurate results. 







On the other side, the accuracy of the experimental drag polars of figure 11 must also 
be questioned. The data were obtained from a two-dimensional insert in a larger wind 
tunnel with no blowing or suction of the wall boundary layer applied. Boundary layer 
separation in the corners formed by the flap and the insert ^wa'lls most, likely tpoks place 
at higher flap deflections and higher angles of attack, thereby causing a , deviation from 
the desired two-dimensional flow pattern. .Many the.mb^ tas^dBsa^qpg puMigned in 
reference 12' support this observation. * ; 1 M ■ ? t.'-'M-iM- 


- w 
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Figure 7. — Lift and Pitching Moment Characteristics of GA (W)- 1 With 30% Chord Flap 






j Figure 8. — Lift Curves of GA(W)-1 With 30% Chord Flap 
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0 {Experiment: Wentz, Seeth aram, & Fiscko 
.... . _ Baseline version, modified cove geometry 


Version C 


FT Fixed transition 
S Separation 
LS Laminar short bubble 



Figure 10. — Surface Pressure of GA(Wj-1 With 30% Chord Flap 








IPressure figure 12 for the 20® flap angle and in figure 10 

for the 30° case. TThg solid lines represent the theoretical predictions of Version C which 
reasonably: well approximate the experimental surface pressures in all cases 
Differences are observed in the cove region and near the upper surface trailing edge. 
Mispredictions are expected from an aerodynamic theory that does not attempt to 
properly model the recirculating flow in this region. No explanation is offered for >th| | 
discrepancy between theoretical and experimental pressure distributions near the uppe&J 
surface trailing edge. Figure 12 further shows theoretical results of the baseline version 
of the code obtained by Wentz, Seetharam, and Fiscko (ref; 13). The cove geometry of 
the GA(W)-1 had been modified by these authors to produce a better match oi theory 
land experiment in the cov^ region. At this pointfit should be emphasized that all other 
comparisons shown in t|is document are based on the true* airfoil geometries without 
modifications! I - * * * * fe 


I 

■i # i- P 




4 % 


The computed boundary layer thickness (Version C) on both surfaces of the main i 
component of the GA(W)-1 high lift airfoil is shown in figure 13. The boiindaryV ll^M f 
lh ^ n ®? s _ 6 u ; 011 . upper surface gro w s continuous ly , whereas, Si on the. lower surface i 
.suddenly increases at the entrance to the cove and later nearly vanishes at the lower 
|sgrface trailing edge. This thickening of the boundary layer has been noticedjin iall 
. theoretical results and seems, at least qualitatively, to be a realistic description : of the 
flow pattern m the cove region. The boundary layer at the lower surface trailing edge is 
too thin which results in erroneous initial values for the subsequent calculation an Me# , 
core region. The potential core size is overestimated at the slot exit which in turn shifts 1 
fo e end of the core reg i on too far downstream. ' •••• < * 




Figure 14 through 19 contain data of the GA(W)-1 with a 30% chord trailing edge flap | 
at 40° flap angle. Experimental data' (ref. 14) clearly shows separation on the flap ri&ll.1 
surface for this configuration. Presumably this is the reason theoretical And ■ 
. ? x P e r i men tal lift coefficients disagMeTfig- 14^Tiej)Eba8ure distributions of figure 15 
* verify the separate^ flow copditiomon^the uppgr surface of Uie flap . In tips area, the 1 
theoretical pressure data dp not agree well with measured pressures. as large scale flow * 
^separation is not modeled by the program. The agreement is better along the surface of 1 
the main airfoil, component, but substantial differences are noted near the upper surface 1 
eading edge aqd, as already observed at lower flap settings, in the cove region. The 

reader should note* that urogram Vpthiah 11 ni*Ailii/tha^m/ivn . . ■ 

second half of 
displacement 
^displacement 

»>-auiuci ixuc w uie ttuiun. lmsresuiteain an 
| unrealistic couplin g of upper and lower surface pressure distributions *that^an cWarly ' 


the upper surface of the' vv 


Figure 16 presents a comparison of skin friction coefficients on ;V 

GA(W)-1. The corresponding predictions of the point of separation ont the* flap upper 
surface are shown in figure 17. Theoretically, ^separglion is not well predicted; which is 
probably due to the misprediction of the pressures on the flap upper surfacd 

i... ’ll i - r rr • v .• XV 


-* ln 2 r 2J ements - ha X e 1° made 'by properly modeling large scale separation 
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© - ' Experiment: Wentz, Seetharam, & Fiscko 

Baseli ne version, modified cove geometry 
1 ; Version C t 

FT Fi xed transition 

S Separation 

; LS Laminar short bubble 



Figure 12. — Surface Pressure of GA(W0-1 With 30% Chord Flap 


25 










/ © i Experiment 
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Figure 14. —Lift Curve of GA(W) 1 With 30% Chord Flap 
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O Experiment: Seetharam and Wenz 



T 


Figure 15. — Surface Pressure of GA(W)-1 With 30% Chord Flap 
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Figure 18 again emphasizes the coupling of upper and lower surface data in jthe, 1 
computations of program Version A. As shown earlier in figures 13 and 15, such a- j 
coupling of upper and lower surface data is not present in the theoretical prediction' of i ] 
Version C. ' ,35-1 




Figure 19 demonstrates the inability of the NASA-Lockheed program to simulate 
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A few results for the two-element airfoil, Foster (ref. 15), are shown in figures 20 and 
21. Theoretical and experimental pressure distributions are only in partial agreement 
and, also higher values of the boundary layer displacement thickness on the upper 
surface of the flap are estimated by Versions A and B of the program. Similar results 
were calculated by program Version C. The reasons for this discrepancy are not known. 
Most likely the correspondence between test data and airfoil configurations in Foster’s 
publication (ref. 15) was misinterpreted since test results in the referenced docume nt do 
not clearly define geometric parameters of the tested configurations (figs. 20 and 21). 
Furthermore, the Reynolds number and the Mach number are not stated and must be 
estimated on the basis of the given tunnel speed (table 2). For these reasons, Foster’s 
high.Jift data are of limited yalue/oij an evaluation pf ^h^ performance predictions’ of j 


the computer program. 
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Lift and drag coefficients of Ljungstrdm’s two-element airfoil (ref. 16) are shown in , 
figure 22. Theoretical predictions of program Version B were only successful a t a zer o 

{degree angle of attack, hut did not. converge to' a soMtiSh-^af^iithWSppel'j 
[was made to u se Ver sion 0 for a i^fformaneg analvsisWLjuhgstrSMs.iaai^feyl-.^ 
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The Boeing four-element % high Jif1|;airfoil (fig. 1 and table 3) 'was used as the , main test 
case for multiple airfoils? Itf* consists of a wing section with a leading edge flap and a 
double slotted trailing edge flap. Global airfoil parameters and detailed distributions of 
surface pressures and boundary layer data are available? for comparisons. The data were 
obtained in the Boeing research wind tunnel (BRWT) on a model with 2 feet unextended 
wing chord and 5 feet span. Careful blowing of the wall boundary layers was applied in 
order to achieve a two-dimensional flow pattern across the whole span of the airfoil. 


The lift curve and the drag polar of this airfoil at a Reynolds number of two million, * 
based on the wing reference chord, are given in [figures 23 and 24. The experimental lift j 
coefficients are ^balance; data that are within 1.5% of the lift obtained hy pressure^! 
integration. The profile "drag of the airfoil is the result of wake rake measurements 
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Figure 19. — Boundary Layer Profiles, Comparison of Theory and Experiment for GA(W)-1 
With 30% Chord Flap, a = 2. 7°, S F = 4(f 






Figure 20. — Pressure Distribution on Flap for Foster Two-Element Airfoil 


r ; (a) Boundary Layer Displacement Thickness 
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Figure 22. — Lift and Drag Comparison for Ljungstrom Two-Element Airfoil 














supporting brackets ancl v "pressure”' taps. The maximum spanwise variation of the 
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downstream of the flap support brackets and the lowest downstream of the pressure^ 


All attempts failed when "using the baseline version of the program and program 
Version A' "to obtain a converged solution for this airfoil. Program Version B arrived at 

[co nverged solutions Wf^^S^^aSd^'O^hgle of attack,Jmt,,yp ^ a\ 

~ -N J Ll ~ a- /r* oo\ Tk „ • : j r n 


considerable amount (fig. 23). The improvement of the predictions by Version C is- 
remarkable. But the reader should note that the potential flow solution already 
provides a very good approximation to the lift curve. Program Version C overpredicts 
the lift at high angles of attack, which is expected since the program does not model 
flow ;S%>‘aration. However, the < program accurately indicates the onset of trailing edge 
stall at about 169 angle >of attack, thereby warning the user that above this angle of 
attac&i tfoperateslputside atsfranige of validity 
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The theoretical values of the profile drag of Version G, (figs. 24 and 25), is relatively 
close to the measured profile drag. In judging the quality of the agreement of the two 
types of drag curves, the reader should recall the problems of two-dimensional high lift , 
testing and the uncertainties in applying the Squire and Young formula to. theoretical^ 



Theoretical pitching moment characteristics are compared with experimental data in 
figure 26. The discrepancy of the curves at higher lift values is due to trailing edge stall 
that is not modeled by the program. • 


Figure 27 and table 4 demonstrate the excellent convergence characteristics ofjtbe new 1 
program Version^C.vTabJe 4 lists the iipdates Qjy^'^psijbi^n'lpif We wakp (pmterlinesl 
emanating from^ t^ilir^^^^^e^i^^^irHc^pp^i^ during^e^h.ccxcle^ o%jtgeJ 
pveraH^itiqia procedu^V ^ s •' 4r 

T ■ : &» ??■■■■'- __ a. 

Figures 28 through 33 contain comparisons of theoretical and experimental surface • 
pressures at 8.4° angle of attack. These figures confirm the earlier findings that Version 
C indeed provides the best theoretical results. Differences between the theory of Version 


C and experiment are noted in cove regions and on the second trailing edge flap (fig. 28,1’ 
32, and 33). Again, the need to model the flow in the cove region of airfoils is apparent. ■■ 
Other differences in the surface pressures are due to an insufficient number of surface 
points used in the discretization of the airfoil geometry. Figure 34 contains the number 
of surface points of each airfoil component. Numerical experience in aerodynamics has 
shown that at least 50 to 60 surface points uare nee1le<? fqr an\ddqu4te‘ representation 
| of each airfoil component. The insufficient discretizatibn of Ifcie^Boeiiig^iijh^lift airfoil 
was dictated by the limitations of program Versions^ A and to a' total <k- 165 ^surface * 
points? Program ^Version Caused the same %umber%f^surface%poin® $>• ensure '#feir 
comparison of the lprbgram' vprsior&; but is not^lifnitedi t^ sSchf gft small number of 
geometry points. % ' h K' ^ *■ 4 

•* .» * ■** 4 ♦„ ft '• <*h ‘^i'% #■'.«; ’t'"' % 

Figure 34 demonstrates the capability of Version C to, compute smooth streamlines in? 
the vicinity of airfoil surfaces. . I..- ■> ~ -0 ; '■ ; !l 
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Figure 27. — Convergence Characteristics of Program Versions for Boeing 
Four-Element Airfoil 



Table 4. — Convergence Characteristics of Version C of 
NASA/Lockheed Program 


Boeing Four-Element Airfoil 


| a = 8.4° 
;M oo = 0.16 
R n = 2 x 10 6 


Iteration number 

Lift 

Wake centerline updates 

1 

2.211 

2 

2 

2.158 

1 

3 

2.176 

1 

4 

2.175 

0 

5 

2.173 

0 




| a= 12.45° 

I Mao= 0.16 
R n = 2x 10 6 


Iteration number 

Lift 

Wake centerline updates 

1 

2.748 

1 

2 

2.676 

1 

3 

2.698 

1 

4 

2.701 

1 

5 

2.699 

1 
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Figure 28. — Slat Surface Pressures of Boeing Four-Element Airfoil 
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Figure 29. 
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Figure 30. - Wing Surface Pressures of Boeing Four-Element Airfoil 


O i Experiment 

— — i Version A , 

— — = i Version B 

— Potential flow 

a = 8.4° 
:M m = 0.16 
R n = 2x 10 6 



Figure 31. — Pressure Distribution Comparison for Main Flap of Boeing Four-Element ■ 
High-Lift Airfoil 
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Figur e -.3 5 through 39 compare predictions^ of* boundary layer' parameters! figure ■38,t- 
deserves special attention. "Boundary layer^velocity profiles or^ tlfe upper surface of the 
main component are shown at ^several chordwise ’statioiis.<*'Th'e experimental- velocity 
profiles reveal that very little confluence *of slat wake arid wing boundary layer has 
taken place, and that' an^iriitially^existing weak eonfluent^boundary layer above the j 
wIm has degenerated early irilo r an ordinary t'urbulenti'boundary^layeijf Thisjfeature of 
’the flow field is very well simulated'by Version C.# ,f .jt J 1 # M & & j# 

' ;v .. ' * " * •* •’ f ' » #■ * •/ #./’# / $ r * # #> . * 

LJUNGSTROM’S AIRFOIL* * # » 0 ,# *>' # # ^ ^ ^ ^ , | 

Fi gur e 40 shows lift and drag curves of the Ljungstrom? four-element airfoil (ref. 16). 
Gpnverged^hririreticaksolufidifslwere^rily obtained^at a = 10° and'l'Gf ‘angle of' attack^; 
Obviously, successful theoretical, predictions, for this fcp>e.$of, airfoil jwill require, a models^ 
C0 T e separation- m # 41 * & M jit B m Jt A Jfr J* ,, 




Region I of confluent boundary layer 



Figure 36. — Viscous Flow on Upper Surface of Boeing Four-Element 
High-Lift Airfoil — a = 8.4° , = 0.16, R^ = 2 x icfi 
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CONCLUSIONS 


Several versions of the ,NASA-Lockheed computer program, were evaluated on the basis 
of experimental high lift data that had not been available when the original version of 
the code was published. Based on a relatively short evaluation phase of a few months, 
the following conclusions about the reliability and quality of the theoretical predictions 
are drawn. " ^ *” * 


Program ‘•Versioif&C, developed by The Boeing Company, is the most reliable and j 
accurate version. It produced converged solutions within a few iteration cycles for f*! 
all test cases run, while the other program versions frequently failed to provide 
any result at all. The improvement of the accuracy of the program predictions isfvf 
due to several modifications of the aerodynamic model; mainly due to a better i 
representation of the boundary layer displacement effects, and an improved model I 
of the potential core region. •* j 

All program versions produced the best results in cases where most of the flow is 
attached to the airfoil’s surface. This is consistent with the basic assumption of || 
attached flow, but the range of applicability of the program should be extended by 0 
adding separation models for the cove region and for trailing edge stall. 

The usefulness of the confluent boundary layer method of Goradia and its 1 
modification utilizing Coles’ velocity profile for the purpose of predicting the onset ■] 
of confluent boundary layer separation, has yet to be tested. Optimized ■! 
configurations were chosen for most of the program evaluation with little 

confluence of wakes and boundary layers. ^ * . * j 

• * ' y s & 4 0 0 „ # . ■ «*.• * \ ■ i 

t ~f s? , t 0 , • *■ 'it 

The performance of the program needs to be tested for configurations at off-design fl 

. • & * # ’ • 

conditions. J m 

m ■ -it 

The evaluation of the computer program was hampered by the shortage of reliable 
experimental high lift data. Additional wind tunnel testing of some of the more 
important high lift airfoil configurations would increase the confidence in their ^ 
performance predictions. * ; <* ■ ■ 


Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle, Washington 98124 
December, 1977 
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